function [elements,processInfo] = InitializeSolutionStep(elements,model,processInfo)
ne = model.countElements();
dim = 2;
for i=1:ne
    e = model.getElement(i-1);
    f = e.getFace();
    ue = f.getUHatLocal();
    nn = f.countModes();
    ng = e.getNumGaussPoint();
    for j=1:ng
        xi = e.getGaussPoint(j-1);
        B = zeros(3,dim*nn);
        B = e.computeB(f,B,[xi(1) xi(2)]);
        %compute strain
        epsilon = B*ue;
        %return map implicitly
        elements(i) = implicitReturnMap(j,elements(i),epsilon);
    end
end
end
